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Abstract - The physical properties of hole-doped cuprate high-temperature superconductors are 
heavily influenced by an energy gap known as the pseudogap whose origin remains a mystery 
second only to that of superconductivity itself. A key question is whether the pseudogap closes 
at a temperature T*. The absence of a specific heat anomaly, together with persistent entropy 
losses up to 300K, have long suggested that the pseudogap does not vanish at T*. However, 
amid a growing body of evidence from other techniques pointing to the contrary we revisit this 
question. Here we investigate if, by adding a temperature dependence to the pseudogap energy and 
quasiparticle lifetime in the resonating-valence-bond spin-liquid model of Yang Rice and Zhang, 
we can close the pseudogap quietly in the specific heat. 


The physical properties of hole-doped cuprate high- 
temperature superconductors are strongly influenced, over 
a wide range of temperature and doping, by a depletion in 
the electronic density of states known as the pseudogap. 
In momentum-space it is manifest as a gapping of the large 
hole-like Fermi surface near the antinodal regions of the 
Brillouin zone, at (±7r,0) and (0,±7r), leaving behind un¬ 
gapped “Fermi arcs” [T]. The origin of the pseudogap re¬ 
mains a mystery second only to that of high-temperature 
superconductivity itself, and it is widely hoped that by in¬ 
vestigating the former we might uncover valuable insights 
for understanding the latter. A key question is whether 
the pseudogap closes at a temperature T*. In recent years, 
evidence has been building that suggests that it does. 
These include abrupt changes in the Kerr Effect [2], time 
resolved reflectivity [3] , as well as the direct observation of 
a reconstruction of the antinodal electronic structure by 
angle-resolved photoemission spectroscopy [3B] (ARPES). 
In this work, we aim to reconcile those results with ther¬ 
modynamic measurements, in particular the electronic en¬ 
tropy and specific heat, which have long suggested that the 
pseudogap is temperature independent [5H8]. 

By way of introduction, the electronic entropy is defined 
as S{T) = -2kB f fw{E, T)N{E)dE [9] where N{E) is the 
density of states, and fw{E, T) is a “Eermi window” which 
expands with temperature and is related to the Fermi dis¬ 
tribution function / by / In / -|- (1 — f)ln{l — /). Put sim¬ 
ply, S{T) is a count of the thermally active states. The 
electronic specific heat coefficient is given by the temper¬ 


ature derivative of the entropy, ^{T) = dS{T)/dT. Three 
apparently universal observations have been made from 
high-resolution differential specific heat studies on a va¬ 
riety of hole-doped cuprates [SHSlflO]. These are: i) a 
loss of entropy in low to slightly overdoped samples, that 
persists right up to the highest temperatures measured. 
The entropy decreases at a rate of about 1 ks per doped 
hole; ii) a collapse in the magnitude of the specific heat 
jump. Ay, at Tc below a critical doping of 0.19 holes/Cu; 
and iii) a smooth downturn in the normal-state electronic 
specific heat with no specific heat jump at T*. These 
features were originally modeled by Loram in terms of a 
temperature-independent non-states-conserving V-shaped 
gap, pinned to the Fermi level {Ep) of a flat density of 
states. The gap widens with reducing doping [8]. In con¬ 
trast to the superconducting gap, where the low-energy 
states are pushed just above the gap edge, it is surmised 
in this model that the pseudogap redistributes those states 
to much higher energies. In this scenario T* represents an 
energy scale where thermal fluctuations become compa¬ 
rable in magnitude to the size of the pseudogap, rather 
than a phase transition temperature. If one tries to fill in 
such a pseudogap with temperature, thereby simulating 
expanding Fermi arcs mm, problems arise. Firstly, the 
lost entropy is eventually recovered, contradicting (i). Sec¬ 
ondly, a kink in the entropy appears at T* together with 
a corresponding jump in the heat capacity m , contra¬ 
dicting (iii). And finally, we might expect a double-peak 
structure to appear in the superconducting anomaly near 
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critical doping, where T* is less than Tc, altering the dop¬ 
ing dependence of A-yiTc) compared to (ii). But perhaps 
this just means that this model is incomplete, and if so, 
what are we missing? 

In the following we will investigate the effects of a tight- 
binding density of states, thermal lifetime broadening, and 
the combination of these with a Fermi-surface reconstruc¬ 
tion model for the pseudogap given by the resonating va¬ 
lence bond spin liquid ansatz of Yang, Rice and Zhang 
(YRZ) [13]. Detailed descriptions of the YRZ model have 
been published several times [IMS, but for complete¬ 
ness we briefly list the equations used in this work. In 
the normal state the coherent part of the electron Green’s 
function is given by 

G{k,u;,x) = -(1) 


where = —2t(a:)(cos kx + cos ky) — 4t'(a:) cos kx cos ky — 
2t"{x){cos2kx + cos2ky) — fip{x) is the tight-binding 
energy-momentum dispersion, = —2t{x){coskx + 
cosky) is the nearest-neighbour term, and i?g(k) = 
[i?°(a:)/2](cosfca; — cosky) is the pseudogap. The chemi¬ 
cal potential Hp{x) is chosen according to the Luttinger 
sum rule. The doping-dependent coefficients are given 
by t{x) = gt{x)to + {3/8)gs{x)Jx, t'{x) = gt{x)t'Q and 
t"{x) = gtix)tQ, where gtix) = 2x/(\ -I- x) and gs(x) = 
4/(1 -I- xY are the Gutzwiller factors. The bare parame¬ 
ters t'/to = —0.3, t'V^o = 0.2, J/to = 1/3 and y = 0.338 
are the same as used previously [T3] . Equation [T] can be 
re-written as 


G(k,a;,x) = ^ 

a—± 


gt{x)W^{x) 

w-E^[x) 


( 2 ) 


where the energy-momentum dispersion is reconstructed 
by the pseudogap into upper and lower branches 




that are weighted by 


Ck + C] 


^k =2 


1 ± 


(a + ?^)/2 


0(ek + C£)/2F + i?,^(k) 


+ El{k) (3) 


(4) 


In the superconducting state there are four energy 
branches EEg = -I- A^, where a = ± and 

Ak = [Ao(x)/2](cosfc3, — cosky) is the superconducting 
gap. The density of states (DOS), from which the entropy 
and heat capacity can be calculated, is given by 


N{uj)= ^ gtix)W^[iuZ)^Siuj-E^) + K)^Siuj + E^)] 

a—±.k 

(5) 

where (ug)^ = 0.5(1 -I- E^/Eg) and (ug)^ = 0.5(1 — 
E^/Eg) are the Bogoliubov weights. 


The reason for choosing this model is because it suc¬ 
cessfully describes experimental data from a wide range of 
techniques [16], including the specific heat [miTi]. How¬ 
ever, the previous works did not consider a temperature 
dependent pseudogap term. In fig. [1] we plot the calcu¬ 
lated energy momentum dispersion in the superconduct¬ 
ing state along cuts in the ky direction near the antinodes 
for a;=0.12, both with, and without {Eg=0) the pseudo¬ 
gap. The results reproduce the ARPES-derived disper¬ 
sions measured below and above T* respectively [111], 
providing compelling evidence for the closure of the pseu¬ 
dogap at T*. Key details are reproduced such as the sep¬ 
aration between the momentum of the minimum binding 
energy of the dispersion ko from the Fermi momentum 
kp^ a signature of non-particle-hole symmetric order [?]. 
Moreover we can identify the flat dispersion of the shoul¬ 
der feature observed in ARPES energy dispersion curves 
[3] as belonging to the Bogoliubov dispersion arising from 

the upper YRZ band, — 

Since we wish to understand the effect of adding a tem¬ 
perature dependence to Eg^ from here onwards we fix the 
tight binding coefficients to their values at a;=0.20 and 
neglect the gtix) prefactor in the equation for the den¬ 
sity of states. (Normally the x dependence of these terms, 
which narrow the bands but reduce the magnitude of the 
DOS, would complement rather than counteract the pseu¬ 
dogap.) To fit the experimental entropy data to is set to 
0.225 eV. Beginning for a moment without the pseudogap, 
the defining feature of the tight-binding DOS is the van 
Hove singularity (vHs) located just below Ep for a;=0.20 
(see fig.[2](a)). Assuming a rigid shift of Ep away from the 
vHs with decreasing doping results in a persistent decrease 
in entropy, as shown in fig. [2](d). However at 300 K, the 
rate of decrease is only 0.33 fc^/hole compared with the 
observed 1 fcs/hole [8]. 

Lifetime broadening can also affect the high- 
temperature heat capacity, and hence the entropy, 
by smoothing features in the DOS m- From resistivity 
measurements [20j we infer a linear-in-temperature scat¬ 
tering rate (inverse lifetime) given by F = O.Olto + fiksT, 
with a slope (3 that increases with decreasing doping. 
The most computationally efficient way of incorporating 
this term is by convolving the DOS with the lorentzian 
r/7r[(w — E)'^ + F^]. Figure[2](b) illustrates the thermally 
broadened vHs at 300 K for a;=0.20 and 0.14 with /3=1 
and 2 respectively. The entropy decrease is now larger at 
0.7 kp/iiole (fig.[2](e)), but it is still not enough, especially 
at low temperatures. This necessitates the incorporation 
of a pseudogap. 

In fig. He) we add a pseudogap for a;=0.14 by setting 
Eg=54 meV. Based on the ARPES results we initially as- 

t ume that the pseudogap closes linearly with temperature 
ccording to EgiT) = Eg — 2kBT. The van Hove sin¬ 
gularity and lifetime broadening effects are also included. 
The calculated entropy compares well with experimental 
data for La 2 -a;Sra;Gu 04 [6], shown in fig. HO, with the en- 
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Fig. 1: (Color online) (a) Spectral weight at the Fermi level calculated from the YRZ model for *=0.12. The black curves 
indicate the back of the nodal hole pocket, as well as the Fermi level crossings when Eg=0. Dispersions calculated along the 
vertical momentum-space cuts are shown in plots (b) to (e), both with (blue line) and without (red line) the pseudogap. In the 
fcj, = TT cut, (b), the dashed line shows the Bogoliubov dispersion from the upper YRZ band (Ao=0.12to). These results closely 
reproduce the experimental data of Refs. [3] and [5] 



Fig. 2: (Color online) Density of states for *=0.20 and 0.14 in the case of: (a) a rigid shift of the Fermi level and no 
lifetime broadening; (b) the addition of thermal lifetime broadening terms nksT and 2'KkBT respectively; and (c) a YRZ-like 
reconstruction. The corresponding electronic entropies are shown in plots (d) to (f). Plot (f) includes experimental data for 
La 2 -xSr 3 ;Cu 04 [6] and the calculated curves include the thermal broadening terms. The two fits to the *=0.14 data correspond 
to the two pseudogap temperature dependences of Eg{T). 
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Fig. 3: (Color online) Electronic specific heat for a YRZ-like 
pseudogap that closes as Eg = Eg — 2kBT in the presence of 
lifetime broadening F = O.Olfo + ksT. Note the absence of a 
specihc heat jump at T*. 


tropy decrease approaching the observed 1 fcs/hole. The 
low-temperature fit can be further improved by taking a 
more gradual initial T-dependence given by 


EgiT) = El 


2 — 1 / tanh 


El\n2,\ 

^keT J 


( 6 ) 
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We now turn to the specific heat coefficient, 7 . Fig- 
ure[3]shows 7 ( 1 ") calculated for a 40meV pseudogap which 
closes linearly with temperature in the presence of lifetime 
broadening with /3=1. (This of course does not contain 
the additional experimentally observed contributions from 
fluctuations near Tc.) There is no obvious jump at T*, 
only a slight change in slope. If Eg(T) was rounded near 
T*, due to doping inhomogeneity for example, the specific 
heat would become even smoother there. Finally, in fig. |T] 
we plot the doping dependence of the specific heat jump 
at Tc assuming a parabolic superconducting gap doping 
dependence A(a;) = 0.103to[l — 82.6(a; — 0.16)^], and the 
YRZ pseudogap doping dependence T°(a:) = 3<o(0.2 — x) 
for X < 0.2. Here we take the closure of the pseudogap to 
lie at X = 0.2 in continuity with YRZ, however it has been 
extensively shown that this occurs at slightly lower dop¬ 
ing X = 0.19 m- The pseudogap model reproduces the 
collapse of the specific heat jump as reported for example 
in refs. m and [5]. Note that here we have taken a dop¬ 
ing independent lifetime broadening, /3 = 1. Increasing /3 
with decreasing doping would increase the rate of collapse 
of A 7 (rc). 

To conclude, the absence of a specific heat jump at T*, 
together with persistent losses in entropy at high tempera¬ 
tures, has long been taken as evidence that the pseudogap 
does not close there. Driven by a growing body of evidence 
from other experimental probes pointing to the contrary 
we have explored this question. By adding a linear-in- 
temperature scattering rate to a YRZ-like reconstruction 
model, it is possible to close the pseudogap quietly in the 
specific heat. A similar result is expected for the antiferro¬ 
magnetic Brillouin-zone-folding Fermi surface reconstruc¬ 
tion model [22) . The entropy recovery expected from the 


Fig. 4: (a) Superconducting electronic specific heat jump for 

a;=0.14, and 0.17 to 0.22, in the presence of a YRZ-like pseu¬ 
dogap that closes as Eg = Eg — 2kBT and lifetime broadening 
r = O.Olto + kBT. (b) Doping dependence of Tc (squares), T* 
(triangles), and the specific heat jump at Tc both with (filled 
circles) and without (empty circles) the pseudogap. 

closing gap is offset by scattering-induced broadening of 
the van Hove singularity. This scenario could be tested 
experimentally by searching for an ongoing divergence be¬ 
tween neighbouring entropy curves above T*. 
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